DC current through a superconducting two-barrier system. 
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We analyze the influence of the structure within a SNS junction on the multiple Andreev reso- 
nances in the subgap I-V characteristics. Coherent interference processes and incoherent propagation 
in the normal region are considered. The detailed geometry of the normal region where the voltage 
drops in superconducting contacts can lead to observable effects in the conductance at low voltages. 
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The physics of SNS junctions of microscopic size have attracted a greal deal of attention recently |Q-g]. A variety 
of experiments can be understood by modelling the constriction by a number of one dimensional superconducting 
h . channels in parallel H. Each channel is described in terms of its transmission coefficient. This model assumes that, 
at finite voltages, the potential drop occurs in a normal region close to the barrier. The size of this region is implicitly 
fixed when specifying the boundary conditions for the quasiparticle wavefunctions in the two superconducting leads. 
The matching conditions used so far are equivalent to the assumption that the width of the normal region is much 
smaller than other relevant length scales in the system, such as the coherence length, £o- If this is the case, Anderson's 
J • theorem can be used to justify the lack of variation of the superconducting properties of the leads in the proximity of 
the constriction H|. 

Corrections to the preceding picture are expected when the potential drop, and the effective barrier region are 
of sizes comparable to the scales which determine the superconducting properties. We can expect that, near the 
constriction of size L, the mean free path of the superconductor, /, will not exceed L. The effective coherence length 
is given by £ = \/£fil ~ \fSjyL- According to Anderson's theorem, the relative corrections to the superconducting 
properties due to imperfections of size d are of order 0{L/£). Thus, the small parameter which justifies the existence 
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[--». of an abrupt barrier at the junction is oc y/ ' L/^q. For contact regions L w lOnm and larger, these effects need not be 

l/"~) . negligible for materials such as Al, and should be more relevant for Pb or Nb junctions. 

In this work, we analyze the leading corrections to the abrupt barrier limit by assuming that the barrier region 
has an internal structure. We allow the transmission of the central region to depend on energy. The simplest such 
situation is to consider that the junction is made up of two barriers, which define the central region of the contact, 
and that that transport is ballistic in between the barriers. Recent experiments suggest that such a geometry can be 
C^ ■ manufactured with existing technologies p-ft§. 

The study of microscopic models of Andreev reflections was started in 1982 by Blonder et al pi, who analyzed 
the current which flows through a normal metal interface. They considered the case in which scattering takes place 
only at the interface. The Landauer's approach is generalized to account for Andreev reflection. An electron incident 
from the normal part reaches the interface and its probability to cross it or been reflected is calculated. Electrodes 
are in equilibrium. Thus the probability of each departing process and the occupation of the states is controlled by 
the Fermi function. Adding all possible processes current is calculated. 

In that reference [Q , it was shown that at large voltages the IV curves corresponding to a NS interface are linear 
with a slope equal to the normal state conductance. However they are displaced with respect to the normal state IV 
curve by an amount called excess current. Al low voltages V < A the shape of the IV curve is strongly dependent on 
the transmission of the interface. At large values of the transmission, the current at low voltages is finite p| . 

Later Octavio et al. |12j, proposed a model to analyze transport in superconducting constrictions. They modeled 
the system as a superconducting-normal-superconducting (SNS) system and analyzed the transport through two SN 
systems connected in series. Normal scattering was modeled with 8 barriers at the NS interfaces. 

Their model is semiclassical. All possible scattering processes are added, weighted by its probability. Probabilities 
and not quantum mechanical amplitudes are considered. The interference between scattering processes is neglected. 
With this model, it was possible to explain qualitatively the appearance of subharmonic gap structure in the IV 
characteristics of superconducting weak links. The computation of IV curves in superconducting constrictions was 
analyzed by Arnold 0,[15| , within the formalism of Green's functions. In the last years, several authors have improved 



[16 18 1 these calculations, have calculated the fully coherent quantum mechanical I-V curves using different approaches 
[19 24 1 and have generalized the calculations to the case in which transport takes place by resonant tunneling [£5],£6| . 

In this paper we generalize previous calculations to the case in which there are two barriers in a one-dimensional 
superconducting constriction, SISIS sytem. Another interesting and related system has been recently studied by 
Ingerman et al. p8[ . In that work the current through a SNINS with finite length normal regions, is analyzed. They 
find that in these systems the subgap current is enhanced in comparison to that of superconducting constrictions. 
The effect is most pronounced in low transparency junctions. 

We will consider two cases: in Section II the case in which quantum-coherence is maintained between the barriers 
p7[ is analyzed. The details of the calculation are given in Appendix A. In Section III we study the case in which 
both barriers are added in series, without considering quantum interference between the scattering processes at both 
barriers. We end with the main conclusions. 



II. QUANTUM-COHERENT PROPAGATION BETWEEN THE 

BARRIERS 

In this section, we generalize the previous calculations to the case in which there are two barriers, separated by a 
distance L, assuming that phase-coherence in maintained between the barriers, see Fig. [jl. Except for the existence 
of these two barriers, the superconductor is perfectly clean. 
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FIG. 1. Model-system considered to analyze transport in a superconducting constriction when two barriers are present. 

Each superconductor is described by Bogoliubov de Gennes(BdG) equations and is assumed to be in equilibrium at 
its own chemical potential. Boundary conditions on the wave function and its derivative must be given at the interface. 
If there is a finite potential drop V between two superconductors, the chemical potential of both superconductors 
is not equal, eV = [12 — Hi- Here i = 1,2 label the left and right electrodes. A common reference level must be 
chosen. When both superconductors are referred to the same chemical potential the order parameter acquires a 
time-dependent phase |21|. Thus the time-dependent BdG equations must be solved. The time dependence of the 
relative phase 5(f) — (f>2 — 4>\ is given by the Josephson relation 
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To visualize the processes which contribute to the current and perform the calculations, it is easier to work in 
a scattering formalism scheme, in which each barrier is inside a normal region (see Fig. [j]), as developed in p9[ . 
This normal region, of length d <C £, with £ the coherence length of the superconductor, can be introduced without 
losing generality as the condition d <C £ makes the superconducting properties of the constriction irrelevant. The 
superconductor acts as a source of quasiparticles. These quasiparticles incide at the barrier and can be transmitted or 
reflected. Moreover, Andreev reflection processes can occur at the normal-superconducting (NS) interfaces. Assuming 



that A <C Hi, the Andreev approximation can be used. For simplicity, in the following, we assume that all the 
superconducting regions are ideal BCS ones and that the order parameter in all regions is equal. 

Andreev reflection at the NS interface is given by the amplitude a(e) of suffering an Andreev reflection process. If 
the superconductor is of the BCS type |lij 

1 f e - S5 n(e)(e 2 -A 2 ) 1/2 I e |> A, 
a[6) ~A\ e-i(A 2 -e 2 ) 1/2 | e |< A [Z) 

In the absence of a magnetic field, which breaks time-reversal invariance, each barrier is characterized by a scattering 
matrix. 

( n u \ 

In the neck, between the barriers, electrons propagate without suffering scattering events, but acquire a phase e lkL , 
where k is the electronic momentum. 

As a result of all multiple Andreev reflection processes (MAR) and their interference with the normal scattering 
ones, the wave function in the normal region, in zone I, II (a and b) and III can be written as |p_9|| 
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where Vi and V% are the voltage drops at the two barriers and a™(e) = a(e + meVi + neV-2) with a(e) the Andreev 
reflection amplitude. Region Il-a refers to the superconducting part situated just after the left barrier, while region 
Il-b refers to that one just before the right barrier. 

The current is a time-dependent quantity, which oscillates with all the harmonics of the Josephson frequency 

loj = 2eV/h. 

n 

The time-dependence of the current arises from the time-dependence of the superconducting phase induced by the 
voltage drop. Here we will be only interested in the dc current (Jo component) which is the one experimentally 
measured. 

The scattering matrices relates electron and hole coefficients at regions I, II-a,b and III. Scattering matrices of 
electrons and holes are related by Sh(e) = S*i(—e). As explained in appendix A the matching conditions lead to a set 
of matrix equations between the coefficients in the wave functions, which can be recursively solved. 

There it is shown that if the barriers are a distance L far apart, the problem can be mapped into another one with 
a momentum-dependent single barrier T e ff given by 

TiT 9 

Teff(k) = -J-?- : (6) 
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Here 6 is defined as 



r[r 2 = \/R\R 2 t 
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where r[ = —r*ti/t* is the amplitude of probability of reflexion for those electrons which incide from the right onto 
the first barrier. This effective transmission would be the same in the case of a normal central region and/or normal 
left and right electrodes. The superconducting nature of the electrodes enters in the calculation of the wave functions' 
coefficients. 

Due to the dependence on the momentum, the transmission depends on energy and the equations of previous section 
must be generalized to the case in which transmission and reflection are energy-dependent. A similar situation was 
found in p8fl . We assume that Andreev approximation still holds, what implies that the chemical potential is much 
larger than any other energy scale in the problem. Assuming a linear dispersion relation (ph can be written as 
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with/= £±-. 

The procedure to calculate the coefficients A n and B n and the current is detailed in appendix A. Only those 
coefficients A n — A 7 7 l n S nm and B n = B^ mn are non-zero. 

If / = (barriers are at the same point) the IV curves are the ones corresponding to a superconducting constriction 
with a non-energy dependent transmission of value 
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and no new features appear. New I-V curves appear when / =/= 0. Then, the transmission depends on energy and 
oscillates between the values T m i n and T max given by 
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(10) 
(11) 



If one of the barriers is equal to unity, the current is controlled by the other barrier and we again recover the single 
constriction case. The effect of a finite value of / will be more important for larger difference between T max and T m i n . 
Note that, in particular, if both barriers are equal, equal to T, T max = 1.0 and T min = T 2 /(2 — T) 2 , which decreases 
with decreasing T. Some examples are shown if Figs. to |[ 
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FIG. 2. IV curve for the case of coherent transport in a superconducting constriction with two barriers with parameters 
Ti = 0.9, T 2 = 0.1 (j) = and several values of / 
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FIG. 3. IV curve for the case of coherent transport in a superconducting constriction with two barriers with parameters 
Ti = T2 = 0.9 <f> = and several values of / 




FIG. 4. IV curve for the case of coherent transport in a superconducting constriction with two equal barriers for the case 
/ = 0.1, (f> — and several values of the transmission. 



In Fig. H curves corresponding to T\ = 0.9, T2 =0.1 and <fi — 0.0 are shown. The maximum and minimum value 
corresponding to these transmissions are T max = 0.184 and T m i n = 0.05. These values are relatively similar and a 
very strong effect of a finite / is not expected. In fact, as / increases the conductance is slightly reduced. The fact 
that the conductance is reduced and not increased is because we have started from the maximum value at / = 0. 



Another example is shown in Fig. g, in which case curves for T\ = T2 = 0.9, <fi = 0.0 and several values of / are 
plotted. Maximum and minimum values are T max = 1 and T rn i„ = 0.67, which are again not too different. The main 
effect of a finite value of / is the suppression of the current at zero voltage. This effect is also observed in Fig. 
where several curves in the case in which both barriers are equal, / = 0.1 and = are shown. In all the cases the 
transmission for / = would be equal to unity. The deviation of this curve is more pronounced for small transmission 
of the barrier. For all transmissions current at zero voltage is completely suppressed. The reason of this suppression 
is that even at V = 0, all the energies contribute to the current and the transmission is not equal to unity at all 
energies. 




FIG. 5. IV curve for the case of coherent transport in a superconducting constriction with two barriers with parameters 
Ti = T2 = 0.1 4> — and several values of / 

The effect of a finite value of / is specially strong in figure g as T m in — 0.002. The maximum value is 500 times 
the minimum one. The effect is very strong even for a very small values of /, as in the case / = 0.02. 

III. INCOHERENT PROPAGATION BETWEEN THE BARRIERS 

The calculation and features of the IV curves in the case in which propagation between the barriers is not coherent 
are completely different to those of previous section. From current conservation I\ — I2 — I, where li is the current 
which cross junction i. In the normal state, in the absence of coherence Ii = giVi determines the voltage drop in each 
contact, Vi — <?i/(<7i + g2)V and, as a result the total current is obtained adding the resistances in series. 

To calculate the IV curves in the superconducting state we proceed in the same way. I(V) — I\(T\, V\) = h(T2, V?) - 
Ii(Ti,Vi) is computed according to the calculations for a constriction with a single barrier of transmission Tj [191. 
Then the voltage in each junction is numerically determined from the current conservation equation. Note that SGS 
will not appear at values V — 2A/n, but at values Vi — 2A/n. 
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FIG. 6. IV curves corresponding to transport in superconducting constrictions when there are two barriers in the structure 
and coherence is lost between the barriers. The figure shows several curves for the case in which both barriers are equal. 

Examples of IV curves determined by this method are shown in Figs. HandM. In Fig. f| several curves, corresponding 
to the case in which both barriers are equal, are plotted. Subharmonic gap singularities appear at V — 4eA/n, as the 



voltage in each barrier is equal to V/2. The curves are equivalent to the IV curve corresponding to the transmission 
of the barrier, but for a doubled voltage. Fig. |7] show curves for the case in which both barriers are different. The 
position at which SGS appears depends on the value of both barriers. 
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FIG. 7. IV curves of a superconducting constriction in presence of two different barriers, when there is no phase-coherence 
between the barriers. 



IV. CONCLUSIONS 

We have studied the influence of the structure of the normal region where the voltage drops in superconducting 
contacts. We present detailed expressions to analyze the interference effects which arise due to the Andreev reflections 
at different positions within the contact. 

The observed I-V characteristics at low voltages can differ from the usually considered single abrupt barrier case if 
the size of the contact, L, is comparable to the effective coherence length, \f£$L, where £o is the coherence length in 
the clean limit. The influence of inelastic scattering in the normal region, leading to loss of coherence, has also been 
investigated. 

In particular, we have derived the equations for the case in which the transmission through the constriction depends 
on energy. We have analyzed in detail the case in which the dependence in energy arises from the existence of two 
barriers in the contact, in between which the phase is maintained. The I-V curves differ from the ones obtained with 
only one non energy-dependent barrier with any transmission. The effect is strongest when both barriers have similar, 
and low, transmission coefficients. In this case, the energy-dependent transmission oscillates between two values very 
different in magnitude. Even when at zero energy the transmission is equal to unity, we find that current at zero 
voltage is suppressed. 

In the case in which propagation between both barriers is incoherent, current is derived from the current conservation 
requirement. The main feature of the I-V curves is the appearance of subharmonic structure at voltages different 
from 2A/n. 
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VI. APPENDIX A 



The coefficients in the wave functions are related by the scattering matrices giving a set of equations which allows 
us to obtain them. For the electronic part these equations are 
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and similar equations for the holes. From it the coefficients in region I and III are related by 
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with Sh = S* t . From the existence of a unique source term, J5™q , it can be shown that A™ and S™ will vanish 
except for n — m. Thus, the superindex can be dropped and the problem is equivalent to one in which there is a 
single barrier with transmission given by (o). Due to the dependence on the momentum, the transmission depends 
on energy. 

In the following, the transmission is characterized by a scattering matrix 



r(e) t(e) 

As usual, transmission and reflection coefficients satisfy 

R(e)+T(e) = l 
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where R(e) = ^(e)! 2 and T(e) = |i(e)| 2 . In the two-barrier case, the dependence in energy comes from the dependence 
in momentum. We assume that momentum does not change when the particle cross the barrier, thus the reflexion 
and transmission coefficient depends on the energy of the originally incident particle. However, in the following we 
include both the possibility that transmission and reflexion coefficients depend, as in our case, only on the energy of 
the original quasiparticle, and then we define 



■"n(e) =r(e) 



(20) 
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and the case in which both depend on the energy of the quasiparticle, e + 2neV , which is crossing the barrier and has 
emerged after 2n Andreev reflections. Then 



r n (e) = r(e + 2neV) 
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t n {e) =t(e + 2neV) 
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i?„(e) and T„(e) are analogously defined. 

As before we can find the relations between the coefficients of the wavefunctions. After some algebra the equations 
which relate the A n and B n coefficients are 



^n( e )^n+l ~ ^-(n+l)(~ € ) a 2n+\l2nA n — 
i n( e ) r *i+l( _e ) a 2n+2-B n+ i - *L( n+ i) ( _e )a2n+l)r* (e)*n B n + t*_ („ +1 ) (-e)a 2 „+i JS n o 



(24) 



and 
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In the case t n {e) = t(e) = t(—e), equation ( J25| ) can be simplified to 
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which is simpler than (pq) and reduces to the expresion obtained in p9| , but including explicitly the energy dependence. 
Including the contributions from quasiparticles incident on the constriction from both superconductors (electrons 
from the left superconductor -with momentum k and probability /(e) and holes from the right one -with momentum 
— k and probability (1 — /(e)) the dc-current, as calculated in zone I, is given by 
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This expression was defined in [|19| and takes into account the contribution of all scattering processes of quasiparticles 
with both spins. One must solve ( p5[ ) or (Ujb and obtain the B n coefficients and, using this solution, determine A n 
from (|J) 

It only rests now to discuss how equation (p5h or (pTJ) can be solved. Let us write these equations in the form pj| 



V n B n+ i + U n B n + H n B n -i — F n S n o (29) 

where V n , U n , H n and F n depend on energy and are given by the coefficients of the corresponding equation (Eq) or 
@. Consider n ^ 0. Then |§ 
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In terms of these factors S„ 
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If n > 0, substituting the definition of S n in (|30J) 
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Analogously, for n < 



and 5' n <o is given from 
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which in term of recurrent fractions is written 
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Thus, from (37) and (|40|), the terms 5 n can be easily evaluated. To calculate Bq, we write ([35]) for n — as 

[y 5i + U a + HoS-x] B = F (41) 

Thus 

n £_2 ('421 

^ 5i + (7 + ff 5_i V ; 

and the rest of coefncents are determined from (pi) and (B4I) . 
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